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Abstract 



In the minimal supergravity model (mSUGRA), as the parameter tan/3 
increases, the charged Higgs boson and light bottom squark masses decrease, 
which can potentially increase contributions from tH^, gbj and Zibj loops 
in the decay b — > 57. We update a previous QCD improved 6 — > 57 decay 
calculation to include in addition the effects of gluino and neutralino loops. 
We find that in the mSUGRA model, loops involving charginos also increase, 
and dominate over tW, tH^, gq and Ziq contributions for tan/3 > 5 — 10. 
We find for large values of tan /? ~ 35 that most of the parameter space 
of the mSUGRA model for < is ruled out due to too large a value of 
branching ratio B(b S7). For /x > and large tan/3, most of parameter 
space is allowed, although the regions with the least fine-tuning (low mo and 
mi/2) are ruled out due to too low a value of B{b — > S7). We compare the 
constraints from 6 — > 37 to constraints from the neutralino relic density, and 
to expectations for sparticle discovery at LEP2 and the Fermilab Tevatron pp 
colliders. Finally, we show that non-universal GUT scale soft breaking squark 
mass terms can enhance gluino loop contributions to 6 — > ,37 decay rate even 
if these are diagonal. 
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I. INTRODUCTION 



Models of particle physics including weak scale supersymmetry (SUSY) are amongst the 
most promising candidates for new physics at the TeV scale. Of this class of models, the 
minimal supergravity (mSUGRA) model stands out as providing one of the most economic 
explanations for the diversity of soft supersymmetry breaking terms in the SUSY Lagrangian 
0. In this model f^, supersymmetry is communicated from a hidden sector (whose dynamics 
leads to the breaking of supersymmetry) to the observable sector (consisting of the fields 
of the Minimal Supersymmetric Standard Model or MSSM) via gravitational interactions. 
With the assumption of canonical kinetic terms for scalars in the Lagrangian, this leads 
to a universal mass mo for all scalar particles at some high scale Q, usually taken to be 
Mgut- At Mgut, gaugino masses and trilinear terms are assumed to unify at mi/2 and 
Af), respectively. These parameters, along with the bilinear soft term B, provide boundary 
conditions for the renormalization group evolution of the various soft terms from Mqut to 
Mweak- Requiring in addition radiative electroweak symmetry breaking leaves a rather small 
parameter set 

'^O; "^1/25 tail P and sign{fi) , (1.1) 

from which the entire SUSY particle mass spectrum and mixing parameters may be derived. 

The flavor changing neutral current decay of the bottom quark 6 — 57 is well known to 
be particularly sensitive to new physics effects. New weak scale particles {e.g., a chargino Wi 
and top squark ij) could give loop contributions which would be comparable to the Standard 
Model (SM) tW" loop amplitude. Measurements from the CLEO experiment 0] restrict the 
inclusive B X^j branching ratio to be B{B X,7) = (2.32 ± 0.57 ± 0.35) x 10""^, where 
1 X 10"^ < B{B X,7) < 4.2 X 10"^ at 95% CL. Many analyses have been performed 
which compare theoretical predictions of SUSY models to experimental results. 

In a previous report |^, predictions of the 6 ^ S7 decay rate were made as functions 
of the mSUGRA model parameter space. In this study, a number of QCD improvements 
were incorporated into the calculation which reduced the inherent uncertainty of the b —>■ s'j 
decay rate predictions due to the QCD scale choice from ~ 25% down to ~ 9%. SUSY 
contributions to the 6 — 57 decay amplitude included tW, tH^ and Wiqj loops. Results 
were presented for tan [3 = 2 and 10, and for both signs of /x. For /i < 0, large regions 
of parameter space were excluded, especially for tan/? = 10. For /i > 0, all the parameter 
space scanned was allowed by CLEO data: in fact, for some ranges of parameters, the model 
predicts values of B{b 57) close to the central value measured by the CLEO Collaboration. 

Recently, sparticle mass spectra and sparticle decay branching ratios in the mSUGRA 
model have been reanalysed for large values of the parameter tan/5 [§. In the mSUGRA 
model, the range of tan/? is typically 1.6 ^ tan (3 < 45 — 50, where the lower limit depends 
somewhat on the precise value of rrit. For tan /5 ^ 5 — 10, 6 and r Yukawa couplings become 
non-negligible and can affect the sparticle mass spectrum and decay branching fractions. 
The upper and lower limits on tan f3 are set by a combination of requiring a valid solution to 
radiative electroweak symmetry breaking, and requiring perturbativity of third generation 
Yukawa couplings between the scales M^eafc and Mqut- Some optimization of scale choice 
at which the one-loop effective potential is minimized was found to be needed in Ref. [§] in 
order to gain stable sparticle and Higgs boson mass contributions. This scale optimization 
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effectively includes some portion of two-loop corrections to the effective potential 0. It was 
shown that the mass of the pseudoscalar Higgs boson A^, and the related masses of and 
H~^, suffer a sharp decrease as tan/3 increases. In addition, the masses of the lighter tau 
slepton fi and bottom squark bi also decrease, although less radically. Naively, one might 
expect corresponding increases in the loop contributions to b —>■ decay involving bi and 
H+. 

Indeed, Borzumati has shown in Ref. [TDI that as mH+ decreases, the charged Higgs 



contribution to 6 57 decay does increase. However, for large values of tan P, the chargino 
loop contributions increase even more dramatically, and dominate the decay amplitude. She 
further notes that at intermediate to large tan (3 values, there is a non-negligible contribution 
from gq loops. 

In this paper, we re-examine constraints on the mSUGRA model from b ^ s'j decay at 
large tan /3. In doing so, we incorporate several improvements over previous analyses. 

• We present our analysis using updated mSUGRA mass predictions for large tan /3, 
using a renormalization group improved one-loop effective potential with optimized 
scale choice Q = ^m^^m^^. We use an updated value of top mass rrit = 175 GeV. 

• We include in this analysis contributions from gqj and ZiQj loops. These contribu- 
tions require knowledge of the full squark mixing matrices, and hence an improved 
calculation of renormalization group evolution of soft SUSY breaking parameters. 

• As in Ref. 0, we include the dominant next-to- leading order (NLO) virtual and 
bremsstrahlung corrections to the operators mediating 6 — > 57 decay at scale Q ~ rrib. 
In addition, we include NLO RG evolution of Wilson coefficients between scales 
and nib. We also include appropriate renormalization group evolution of Wilson coef- 
ficients at high scales Q > Mw for tW, tH^ and Wiqj loops following the procedure 



of Anlauf [TT[. The corresponding RG evolution of Wilson coefficients for gqj and Z^qj 



loops is not yet available. 

• We compare our results to recent calculations at large tan/? of the neutralino relic 
density and direct dark matter detection rates for the mSUGRA model. 

In Sec. II of this paper, we present some details of our calculations, especially those 
regarding the inclusion of gqj and Ziqj loops. In Sec. HI, we present QCD improved results 
for the b ^ s'y branching fraction in mSUGRA parameter space for large tan/5. In this 
section, we also make comparisons with cosmological and collider search expectations. In 
Sec. IV, we relax some of the assumptions of the mSUGRA framework to see whether gqj 
loops can become large or even dominant. This question is important when considering the 
model dependence of our results. 



II. CALCULATIONAL DETAILS 

The calculation of the width for b —>■ S'j decay proceeds by calculating the loop interaction 
for b ^ sj within a given model framework, e.g., mSUGRA, at some high mass scale 
Q ~ Mw, and then matching to an effective theory Hamiltonian given by 
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H^ff = -^v,,v*Y.c^iQmQ), (2.1) 

where the Ci{Q) are Wilson coefficients evaluated at scale Q, and the Oi are a complete 



set of operators relevant for the process b 57, given, for example, in Ref. |12]. All 
orders approximate QCD corrections are included via renormalization group resummation 
of leading logs (LL) which arise due to a disparity between the scale at which new physics 
enters the 6 — > 57 loop corrections (usually taken to be Q ~ M^), and the scale at which 
the b ^ S'-y decay rate is evaluated {Q ~ ruf,). Resummation then occurs when we solve the 
renormalization group equations (RGE's) for the Wilson coefficients 

Q-^Q{Q)=^,,C,{Q), (2.2) 

where 7 is the 8x8 anomalous dimension matrix (ADM), and 

The matrix elements of the operators Oi are finally calculated at a scale Q ~ and multi- 
phed by the appropriately evolved Wilson coefficients to gain the final decay amplitude. The 
dominant uncertainty in this leading-log theoretical calculation arises from an uncertainty 
in the scale choice Q at which effective theory decay matrix elements are evaluated. Varying 
Q between ^ to 2mb leads to a theoretical uncertainty of ~ 25%. 

Recently, next-to-leading order QCD corrections have been completed for 6 — > 57 decay. 



These include i) complete virtual corrections [|l3l to the relevant operators O2, O7 and O, 



which, when combined with bremsstrahlung corrections [T^JT^ results in cancellation of 
associated soft and coUinear singularities; ii) calculation of 0{al) contributions to the ADM 



elements 'y^j^ for i,j = 1 — 6 (by Ciuchini et al. |]T5|), for i, j = 7, 8 by Misiak and Miinz [jT6 



and for 737'' by Chetyrkin, Misiak and Miinz fT^. In addition, if two significantly different 
masses contribute to the loop amplitude, then there can already exist significant corrections 
to the Wilson coefficients at scale M^y. In this case, the procedure is to create a tower of 
effective theories with which to correctly implement the RG running between the multiple 
scales involved in the problem. The relevant operator bases, Wilson coefficients and RGE's 
are given by Cho and Grinstein for the SM and by Anlauf [1T| for the MSSM. The latter 



analysis includes contributions from just the tW ^ tH and tiWj loops (which are the most 
important ones). We include the above set of QCD improvements (with the exception of 
727"*, which has been shown to be small [O) into our calculations of the 6 — S7 decay rate 



for the mSUGRA model. 

The contributions to Ci{M]y) and Cgi^M^y) from gq and Ziq loops (SUSY contributions to 
C2{M]y) are suppressed by additional factors of g"^) have been presented in Ref. 0, although 
some defining conventions must be matched between Ref. and Ref. [|l^ and Ref. [[13 



The only complication is that the squark mixing matrix T which enters the couplings must 
be derived. To accomplish this, we incorporate the following procedure into our program 
for renormalization group running. 

• We ffist calculate the values of all running fermion masses in the SM at the mass scale 
Mz- From these, we derive the corresponding Yukawa couplings h^, and h^. for 
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each generation, and construct the corresponding Yukawa matrices {hu)ij, {hd)ij and 
{he)ij, where i,j = 1,2,3 runs over the 3 generations. We choose a basis that yields 
flavor diagonal matrices for {hi^jij and {hf,)ij, whereas the CKM mixing matrix creates 
a non-diagonal matrix |[19[| . 

• The three Yukawa matrices are evolved within the MSSM from Q = Mz up to Q = 
Mqut and the values are stored. We use 1-loop RGEs except for the evolution of 
gauge couplings. 

• At Q = Mgut, the matrices {Ahu)ij, {Ahdjij and {Ahe)ij are constructed (assuming 
A{Mgut) = ^0 X !)• The squark and slepton mass squared matrices {M^)ij are also 
constructed, where k = Q,u,d,L and e. These matrices are assumed diagonal at 
Q = Mgut with entries m^dij. 

• The {Ah)ij and {M^)ij matrices are evolved along with the rest of the gauge/Yukawa 
couplings and soft SUSY breaking terms between Mz and Mgut iteratively via Runge- 
Kutta method until a stable solution is found. The entire solution requires the simul- 
taneous solution of 134 coupled RGE's (with some slight redundancy). We use 1-loop 
RGEs except for the evolution of gauge couplings. 

• At Q = Mz, the 6x6 cZ-squark mass squared matrix is constructed. Numerical 
diagonalization of this matrix yields the squark mass mixing matrix F which is needed 
for computation of the gq and Ziq loop contributions. 

At this point, the Wilson coefficients C'r{M]y) and Cs{M\y) can be calculated and evolved 
to Q rrih as described above, so that the b ^ sj decay rate can be calculated 0|. 

As an example, we show in Fig. |I| the calculated contributions to the Wilson coefficient 
Cji^Mw) versus tan/5 for the mSUGRA point mo,mi/2 = 100,200 GeV, Aq = and /i > 0. 
In frame a), we show contributions from WiQ loops, as well as from tW and tH~ . The 
tW contribution is of course constant, while the tH~ contribution is of the same sign, and 
increasing slightly in magnitude. The various contributions from chargino loops increase 
roughly linearly with tan (3 at a. much faster rate, and thus form the dominant components 
of the b ^ sj decay amplitude. In the case shown, there are several large negative as well as 
positive contributions, so that significant cancellations take place. The sum of all chargino 
loop contributions is shown by the dotted curve. In frame b) we show the contributions 
to Ci{Mw) from different gq loops. These contributions vary with tan/? as well and are 
comparable to corresponding contributions from chargino loops. The sum of all gluino loop 
contributions is shown by the dotted curve; in this case, however, the cancellation amongst 
the various loop contributions is nearly complete. 

In frame c), we show the individual and summed contributions from Z^q loops. These also 
increase with tan f3 but, as expected, are tiny compared to the Wiq and gq loop contributions. 
The sum is again shown by the dotted curve. Here, the cancellations are not as complete 
as in the gluino loop case due to the Higgsino interactions of the neutralinos which increase 
with tan /?. 

III. NUMERICAL RESULTS AND IMPLICATIONS 
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A. Constraints from b ^ s-y decay 

Our first numerical results for B{b 57) decay are sliown in Fig. wliere we plot tlie 
brandling fraction versus tan/? for the mSUGRA point mo,mi/2 = 100,200 GeV, Aq = 
and for a) fi < and b) fi > 0. The SM value, after QCD corrections, is B{b —>■ 57) = 
3.2 ± 0.3 X 10~^, where the error comes from varying the scale choice ^ < Q < 2m{, |]7|,p0|. 
The SM result is denoted by the dot-dashed line, and of course does not vary with tan /3. If 
we include in addition the contribution from the tH~ loop, then we obtain the dotted curves, 
which always increase the value of B{b 57). For this parameter space point, including the 
tH~ loop always places the value of B{b S7) above the CLEO 95% CL excluded region 
of B{b 57) < 4.2 X 10"^ 

If we include the full contribution of supersymmetric particles to the computation of 
B{b — > S7), then we arrive at the solid curves in Fig. ^. For the < case in frame a), 
the SUSY loops increase the branching fraction, which increases with tan/3, so that the 
CLEO restriction on B{b — > 57) severely constrains the mSUGRA model for large tan/3. 
For the frame b) case with /i > 0, the SUSY loop contributions generally act to decrease 
the branching fraction, so that much of the parameter space is allowed for moderate values 
of tan/3. Ultimately, as tan/? increases, the decrease in B{b — >• 57) becomes so severe 
that the mSUGRA model becomes in conflict with the CLEO lower 95% CL bound that 
B{b 57) > 1 X 10^^, so that for this particular mSUGRA point, all values of tan/3 > 21 
are excluded for the particular choice of mSUGRA parameters. 

In Figure ^, we show the main result of this paper: the contours of constant B{b 57) 
in the ttiq vs. mi/2 parameter plane for large tan/3 = 35, for Aq = and for a) < and 
b) /i > 0. The contours are evaluated at a renormalization scale choice Q = mi,. The region 
marked by TH is disallowed by theoretical considerations: either electroweak symmetry is 
not properly broken (the large mo, small mi/2 region) or the lightest neutralino Zi is not the 
lightest SUSY particle (LSP). For small mo, the light tau slepton fibecomes so light that 
in the TH region, m^^ < m-^^. The region denoted by EX is excluded by LEP2 constraints 
which require that the light chargino mass fn^^ > 80 GeV. 

In frame a), we see the value of B{b S7) is large throughout the entire parameter 
space plane. The region with small values of m^ and mi/2 which is most favored by fine- 
tuning considerations is in the most severe violation of the CLEO constraint. The 
region below the dotted contour is in violation of the CLEO 95% CL bound for all choices of 
renormalization scale ^ < Q < 2mi,. Thus, SUSY models allowed by CLEO for tan/3 = 35 
and /i < would be required to have mg > 1470 GeV (at mo,mi/2 = 1000,600 GeV) and 
mg > 1600 GeV (for mo,mi/2 = 320,800 GeV). 

In frame b) for /i > 0, we see that the values of B{b ^ S7) are uniformly below the SM 
value, and so usually in better agreement with the CLEO measured value of 5(6 — >■ 57) = 
2.32 ± 0.67 X lO"'^ (where errors have been combined in quadrature). In fact, we note that 
the region with mi/2 — 500 GeV agrees with the CLEO central value for B{b ^ 57)! This 
region corresponds to parameter space points with mg ^ 1200 GeV, and — 400 GeV. 
In this frame, the entire plane shown, except the region below the dotted contour, is allowed 
by the CLEO constraint. The region below the dotted contour falls below the CLEO 95% 
CL value of B{b 57) > 1 x 10~^ for all values of scale choice ^ < Q < 2mb. This is again 
the region most favored by fine-tuning. In this plane, mg ^ 525 GeV and mg ^ 575 GeV. 
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Up to this point, we have only shown results for a constant value of Aq = 0. In Figure 
we show contours of constant B{b — >• 57) in the tjiq vs. Aq plane for mi/2 = 200 GeV, for 
tan/? = 35 and a) fi < and b) fi > 0. For frame a), we see that the branching fraction 
can change by typically a factor of 2 over the parameter range shown, with most of the 
variation occuring for changes in mo, instead of with Aq. The entire plane shown in frame 
a) is excluded by the CLEO bound. In frame b), for yU > 0, the branching fraction can 
change by up to a factor of ~ 4 over the plane shown, again with most of the variation 
coming due to changes in rriQ. In this case, the region to the left of the dotted contour is 
excluded by the CLEO bound for all choices of renormalization parameter ^ < Q < 2mb. 



B. Comparison with relic density and direct detection rates for neutralino dark 

matter 



An important constraint on the mSUGRA model comes from implications for the relic 
density of dark matter in the universe. The idea here is that in the very early universe, 
the LSP (the lightest neutralino) was a constituent of the matter and radiation assumed to 
be in thermal equilibrium at some very high temperature. As the universe expanded and 
cooled, the LSP's could no longer be produced, although they could still annihilate with 
one another. Upon further expansion, the neutralino flux dropped to such low levels that 
further annihilations would rarely occur, and a relic abundance of neutralinos was locked in. 
These relic LSP's could make up the bulk of dark matter in the universe today. 

The neutralino relic density is calculable as a function of mSUGRA model parameter 
space [^. The relic density is usually parametrized in terms of Qh"^, where Q = p/pc, p is 
the relic density, pc is the critical closure density of the universe (p^ = and H = lOOh 

km/sec/Mpc is the scaled Hubble constant with 0.5 ^ h ^ 0.8. A value of Qh"^ > 1 implies 
a universe with age less than 10 billion years, in conflict with the ages of the oldest stars. If 
Qh"^ < 0.025, then the relic density of neutralinos cannot even account for the dark matter 
required by galactic rotation curves. Some popular cosmological models that account for the 
COBE cosmic microwave background measurements as well as structure formation in the 
universe actually prefer a mixed dark matter (MDM) universe with a matter density 
ratio of 0.3/0.6/0.1 for a hot dark matter/cold dark matter /baryonic matter mix. In this 
case, values of flh'^ ~ 0.15 — 0.4 are preferred. 

In Fig. ^, we show contours of constant relic density Qh"^. The region to the right of the 
solid Qh"^ = 1 contour is excluded by the age of the universe constraint, while the region 
below the dotted contour has Qh"^ < 0.025. The region betwen the dashed-dotted contours is 
favored by a MDM universe. The region excluded by CLEO data is below the solid contour 
labelled b — >■ 57. In frame a), we see that combining the two constraints allows only a small 
patch of allowed parameter space with mi/2 > 700 GeV and 300 ^ mo ^ 600 GeV. Over 
almost all of this region, the relic density Qh"^ > 0.4. In frame b), the b —>■ sj excluded 
region hardly intersects with the MDM region, so that large regions of parameter space are 
favorable for cosmology as well as for CLEO constraints! 

In Fig. 1^ we also plot one contour for expected rates for direct detection of neutralino 
dark matter via cryogenic dark matter detectors . The calculations have been performed 
for neutralino scattering from a "^^Ge detector. Current experiments are sensitive to detec- 
tion rates of 1-10 events/day/kg of detector. The goal of such experiments is to achieve 
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a sensitivity of ~ 0.01 events/kg/day by about the year 2000. Towards this end, we show 
the 0.01 event/kg/day contour in the mo vs. mi/2 plane; below the contour the event rates 
exceed the 0.01/kg/day benchmark. In frame a), we see that the region accessible to direct 
neutralino detection coincides with the region with very large B{b s'y) rates well beyond 
the CLEO 95% CL limit. However, in frame b), for /i > 0, there exists a significant region 
with large direct detection rates, which is allowed by the CLEO constraint, and is also in 
the favorable cosmological region! 

C. Implications for collider experiments 

The LEP2 e^e~ collider is expected to reach a peak CM energy of ~ 195 GeV, which 
should allow SM Higgs bosons of mass thhsm ~ 9^ GeV to be explored. For the tan /5 = 35 
value shown in Fig. |^, the light Higgs scalar nih ^110 GeV over the entire mo vs. mi/2 
plane shown. Hence, for these values of tan /5, we would expect no Higgs signals to be seen 
at LEP2. The reach of LEP2 via fifi and WiWi searches is shown by the dashed contour 
just above the region marked EX. This contour is defined by requiring m^^ = 95 GeV and 
m^j = 85 GeV. For both cases of /i < and /i > shown in Fig. |^, the LEP2 sparticle 
reach falls below both the B{b 57) excluded region, and below the Qh'^ = 0.025 contours. 
If we accept the mSUGRA model literally, then the prediction is that LEP2 should see no 
evidence for either a Higgs or SUSY if tan (3 ~ 35. 

The Fermilab Tevatron pp collider is expected to operate at -y/i = 2 TeV in Run 2, and 
to amass ~ 2 fb~^ of integrated luminosity by use of the Main Injector (MI). Ultimately, 
experiments hope to acquire ~ 25fb^^ of integrated luminosity under the TeV33 program. 
Recently completed calculations of the reach of the Tevatron MI for mSUGRA at large 
tan/3 ~ 35 show a maximal reach in mi/2 to ~ 150 GeV in the ^j^+jets channel |^. A 
similar reach has been calculated for TeV33, and finds points with mi/2 — 175 accessible. 
Comparing these regions to Fig. ^ shows that, like LEP2, the reach of Tevatron MI and 
TeV33 are below both the B(b ^ 57) excluded contour and the VLh'^ = 0.025 contour, 
making discovery of SUSY particles highly unlikely for mSUGRA if tan/5 is large. Over 
much of the parameter space plane in frame b) of Fig. |^, however, m/j ^120 GeV, which 
(optimistically) corresponds to the maximal reach for h at TeV33. Hence, if mSUGRA is 
correct and tan f3 ~ 35, then TeV33 experiments may see a hint of the Higgs boson in their 
data sample. Of course, the entire large tan f3 parameter space shown should be easily visible 
in at least the jets+-^j, channel at the CERN LHC, even with modest integrated luminosity. 

IV. MODEL DEPENDENCE OF gq LOOP CONTRIBUTIONS 

It is well known that within the mSUGRA framework the chargino loop gives the dom- 
inant SUSY contribution to the amplitude for the decay b —>■ sj. This can also be seen 
from Fig. |l] where we see that while the contributions from the gluino loops are individu- 
ally comparable (or even larger!) than those from chargino loops, these cancel out almost 
completely leaving only a small residual contribution. In contrast, while there is indeed 
considerable cancellation amongst the various chargino contributions, there is nonetheless a 
sizeable residue that remains. We may understand the large cancellations among the gluino 
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contributions in analogy with the famihar GIM cancellation in the SM: indeed such a can- 
cellation would be exact if squarks were precisely degenerate {i.e., the squark mass matrix is 
proportional to the unit matrix) because we can then, by a unitary transformation, align the 
squark and quark mass matrices, so that the gluino-squark-quark vertex is exactly flavour di- 
agonal. Within the mSUGRA framework with universal soft breaking squark mass matrices 
at the unification scale, squarks are indeed (approximately) degenerate, and gluino loop con- 
tributions to the flavour violating 6 — > 57 decay are suppressed. The GIM-like cancellation 
that we have described above does not occur when Yukawa couplings enter the calculation 
as occurs, for instance, via the higgsino components of chargino and neutralino loops. In 
this case, as can be seen from Fig. |I|, the cancellation is indeed incomplete (particularly for 
the chargino where the large top quark Yukawa coupling enters). 

These considerations lead us to examine whether the breaking of the degeneracy of the 
soft SUSY-breaking squark masses at the unification scale so strongly upsets the delicate 
cancellations that it results in large gluino contributions to the amplitude for 6 — >■ 57 decay. 
Of course, by allowing soft-breaking mass squared matrices with arbitrary off-diagonal en- 
tries, it should be possible to get very large flavour violating gluino interactions. The issue 
that we address, however is whether large gluino contributions are possible even if we choose 
these soft squark matrices to be diagonal at the unification scale. As we will soon see, the 
physics of this ansatz is basis-dependent. 

To parametrize the breaking of the squark degeneracy we begin by noting that we may 
always choose a quark basis so that either the down or the up type Yukawa couplings are 
diagonal at the weak scale. We will call these the d- and u- cases, respectively. Next, these 
couplings are evolved to Mqut, where both up and down Yukawa matrices have off-diagonal 
components. In the d-case, the down type Yukawa matrices get off-diagonal contributions 
just from the RGE, while the up type Yukawas start off off-diagonal right at the weak scale; 
in the u-case, the situation is reversed. Up to now, the choice to work in the u or d cases 
is purely a matter of convention, and indeed in previous sections we have used the d case. 
This is, however, no longer the case if we further assume that the soft breaking squark mass 
squared matrix is diagonal (but not a multiple of the identity) at the GUT scale. This is 
because the transformation that takes us from the ci-case to the w-case does not leave the 
squark mass squared matrix diagonal (except in the case when this matrix is mg x 1). The 
u- and d- cases are thus physically distinct. 

We have, therefore, studied these two cases separately. To keep things simple, we split 
only the 6-squarks (t^ splits with bi, of course) keeping the others degenerate at mg. The 
splitting is given by a single parameter 

X = (m^/mo)^ 

where and mo are soft squark masses at the GUT scale. Thus x = 1 corresponds to the 
mSUGRA case. We further consider three possibilities where (i) just bi, (H) just b^ and 
{Hi) both bi and b^ masses are split from those of other squarks. 

The results of our calculation of the gluino contribution to Cii^Mw) where the squark 
degeneracy is broken as described above is shown in Fig. ^ for the d-case labelled d-diagonal, 
and for the u-case, labelled w-diagonal. In each frame, we have three curves labelled L, R 
and LR for the cases where just left, just right, and both left and right sbottom soft masses 
are different from mo. In our calculation, we have chosen mo = 500 GeV, mi/2 = 200 GeV, 



9 



tan /3 = 35, = and /i > 0. We choose a large value of mo so that the squark masses are 
not dominated by rrig (in which case splitting due to non-universal soft mass term would be 
unimportant). Also, for the large value of tan/5 the bottom Yukawa coupling is significant. 
The following features are worth noting. 

• For X = 1 the value of Cf^^^° {My/) is the same for the up and down cases for reasons 
that we have already explained. 

• For non-degenerate squarks, gluino loop effects are significantly larger in the w-case. 
This may be understood if we recall that in the d-case, the mixing of down type 
Yukawas at Mqut arises only due to RGE. 

• Non-degeneracy effects when just the right squarks are split show only small variation 
with the non- degeneracy parameter x because flavour mixing in the right squark sector 
is suppressed. 

• Somewhat surprising is the fact that despite the large degree of non-degeneracy, the 
gluino contribution to Cj^Mw), which increases by up to an order of magnitude relative 
to that in the mSUGRA case, never becomes really large. For our choice of parameters, 
this may be partly due to the fact that gluinos and squarks are significantly heavier 
than W and t. This is not the complete reason though. In our computation we 
find that the three main contributions (from the two sbottoms and the sl) cancel 
one another leaving a remainder that is typically smaller than 10-15% of the largest 
contribution. While this cancellation is much less complete than in the mSUGRA case, 
we are unable to give a simple argument for why such a cancellation occurs. 

Finally, we comment on the neutralino contributions for the case of non-degenerate 
squarks. We have already noted that cancellations among various chargino contributions 
are incomplete because of the effect of large higgsino couplings to the tt system. For large 
values of tan /?, we may expect a similar effect for neutralinos. We may further guess that 
this effect is largest in the w-case where flavor mixing does not originate solely in the RGE. 
We have checked that for values of parameters in Fig. ^ above a small value of x = 0.0625, 
the neutrahno contribution to Cj{M\y) is indeed enhanced by a factor of 6-7 above its 
mSUGRA value, and further, that this enhancement is largely due to incompleteness in the 
cancellation between various contributions. For very large values of a;, C-j{Mw) is about 1.5 
times its mSUGRA value, but opposite in sign. We thus conclude that while the neutralino 
contribution is somewhat sensitive to the splitting of squark masses, it never appears to 
become very dominant. 

To summarize the results of this Section, we see that with our assumptions, gluino 
contributions to the amplitude for the 6 — > 57 decay never dominate the SUSY contribution. 
This contribution may nonetheless be non-negligible even if gluinos and squarks are well 
beyond the reach of the Tevatron (and its proposed upgrades) as seen in Fig. |^. We emphasize 
though that our conclusion is special to models where all the flavour violation in the gluino- 
squark-quark vertex at the GUT scale comes from non-diagonal Yukawa interactions. Larger 
contributions from gluino loops may be possible in other models. 

Note added: After completion of this manuscript, a related paper by Blazek and Raby 
appeared on the topic of 6 — > 37 constraints on 5*0(10) SUSY models at large tan/? p6| . 
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Since Ref. [|2^ adopts a particular 5*0 (10) framework and does not include the radiative 
electroweak symmetry breaking constraint, comparison of results between the two papers is 
not straightforward. Also, we have not included — cl mixing included in Ref. 
evaluation of the chargino loop. 
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FIGURES 



FIG. 1. We plot the value of the Wilson coefficient C7(M^) versus tan/?, where 
mo,mi/2 = 100,200 GeV, Aq = and /j. > 0. In a), we show the contribution from the tW 
and tH^ loops, as well as contributions from various loops containing charginos. The q contribu- 
tion includes both cl and ul squarks. In b), we show the corresponding contributions from loops 
containing gluinos, and in c) we show the contributions from loops including neutralinos. 

FIG. 2. We show the branching ratio B{b s^) versus tan/3 for the same parameter space 
point as in Fig. 1, but with a) /j, < and 6) /x > 0. The curves shown are for Standard Model pre- 
diction (dot-dashed), SM plus charged Higgs prediction (dots) and the complete SUSY calculation 
(solid). 

FIG. 3. Plot of contours of constant branching ratio B(b sj) in the mg vs. plane, 
where tan/3 = 35, = and mt = 175 GeV. Each contour should be multiplied by 10^^. Frame o) 
is for /X < and frame b) is for /x > 0. The regions labelled by TH (EX) are excluded by theoretical 
(experimental) considerations. The EX region corresponds to the LEP2 limit of > 80 GeV 
for a gaugino-like chargino. 

FIG. 4. Plot of contours of constant branching ratio B{b sj) in the rriQ vs. Aq plane, where 
tan/3 = 35, mi/2 = 200 GeV and nit = 175 GeV. Each contour should be multiplied by 10~^. 
Frame a) is for /x < and b) is for /x > 0. 

FIG. 5. Plot of contours of constant neutralino relic density in the mo vs. mi/2 plane, 
where tan/3 = 35, = and mt = 175 GeV. We also show the region excluded by CLEG data 
on B{b 57) searches (the region below the b ^ sj solid contour). The region accessible to 
direct neutralino dark matter detectors is below the 10~^ contours. The region accessible to LEP2 
sparticle searches is below the dashed contour in the lower- left. 

FIG. 6. Contributions to C7{M\y) from the gluino- squark loop for the d-case and u-case 

2 

mr 

discussed in the text as a function of the non-degeneracy parameter x = —j. The curves labelled 
L, R and LR refer to the cases where just m^^ , m^^ and both mj^^ ^ are different from the universal 
squark mass mo at the unification scale. For reference, we remind the reader that Cf^ = —0.23 
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